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We derive hierarchies of separability criteria that identify the different degrees of entanglement 
ranging from bipartite to genuine multi-partite in mixed quantum states of arbitrary size. 
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Quantum coherence is deemed responsible for a large 
variety of features, ranging from fundamental physical ef- 
fects such as super-fluidity, via a broad range of counter- 
intuitive interference and correlation phenomena with 
potential implications in the realm of quantum informa- 
tion technologies [1] to transport processes [2] even at the 
mesoscopic scale on the border between physics, chem- 
istry and biology [3]. A central coherence property of 
e.g. a photonic wave packet is its coherence length, but 
the extension of such a concept to composite quantum 
systems is by no means straight forward. 

Entanglement theory promises an accurate character- 
ization of coherence properties in multi-partite quantum 
systems in terms of fc-partite entanglement, i.e. the min- 
imum number k of entangled components necessary to 
describe an n-partite system (in literature also referred to 
as depth of entanglement [4] or fc-producibility [5]). The 
definition of these concepts (as given in Eq. (1) below) 
is rather elementary, but, due to its non- constructive na- 
ture, the identification of fc-partite entanglement in given 
mixed quantum states is a largely open problem; up to 
now the theory of bipartite entanglement {i.e. k — 2) 
has been developed fairly well [6], and there has been 
substantial progress in the identification of genuine n- 
partite entanglement {i.e. k = n) [7-9]. On the scales 
in-between, for n > k > 2, however, only punctual knowl- 
edge, typically for states of specific type or size, is cur- 
rently available [4, 5, 10]. 

The ability to probe these scales in-between is highly 
desirable for various reasons: while it is well established 
that quantum computations with pure states necessarily 
require a large amount of entanglement in order to per- 
form beyond the classically achievable [11], the situation 
is not as evident for mixed states as they would occur 
in realistic implementations, since also mixed separable 
states can lead to improved computational power [12, 13]. 
The possibility to identify entanglement properties in a 
more fine-grained version than currently possible for the 
mixed case would certainly help to understand which spe- 
cific features of multi-partite quantum states are really 
necessary for the appraised quantum speed-up. 

In precision interferometry, the full enhancement of 
precision based on n particles can be obtained only for 
a genuinely n-partite entangled state [14]. Entanglement 
between fewer components will result in a precision closer 
to the achievable with n independent particles: identi- 
fying the largest fc-partite entanglement (for k<n) that 



can be realised at given experimental conditions provides 
therefore very rigorous limitations to the achievable pre- 
cision. 

Similarly, such an assessment permits to estimate the 
number of nodes over which coherence in a computa- 
tional network has been achieved [15]. Fast excitation 
transport through molecular- or spin- networks has been 
shown to be associated with quantum coherence between 
an intermediate number of nodes [16], and such coher- 
ence can be identified through fc-partite entanglement 
after projection onto the single-excitation subspace [17]. 
This provides a very accurate characterisation of the spa- 
tial extent over which a multi-partite system displays 
quantum mechanical features, and the environmentally 
induced degradation of coherence can then be followed 
to monitor the emergence of classicality in a rather de- 
tailed fashion. 

Our goal in the present contribution is, therefore, to 
provide for any system size n a full hierarchy of separa- 
bility criteria to characterize multi-partite entanglement: 
the criterion at the top of each hierarchy identifies gen- 
uine rt-partite entanglement, followed by criteria that are 
positive only for states with at least fc-partite entangle- 
ment for fc ranging from n — 1 to 2. 

Before introducing our framework, let us review briefly 
the necessary formal background. A pure state |^'n,n) 
of an n-partite quantum system is considered n-partite 
entangled if there is no separation of the sub-systems 
into two groups, such that \^n,n) could be described as 
the tensor product of states of these two groups. Anal- 
ogously, an n-partite state |4'fe,„) is considered fc-partite 
entangled if it can not be described without an at least 
fc-partite entangled contribution. If a state is not at least 
bipartite entangled, then it is separable. For pure states, 
definition and identification of fc-partite entanglement is 
rather straight forward, but the situation changes dras- 
tically for mixed states: a mixed n-partite state is con- 
sidered fc-partite entangled if it can not be expressed as 
a statistical mixture 
fc— 1 

of at most (fc— l)-partite entangled states with pij > [5]. 
This leads to a rather intricate structure of multi-partite 
entanglement as sketched in Fig. 1. 

The task of our present hierarchies of separability cri- 
teria is to provide a potentially accurate identification of 
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FIG. 1. Simplified, schematic structure of 4-partite states: 
the grey circles depict quantum states that are separable 
with respect to three different bipartitions (the biseparation 
ABC ~ D and biseparations into pairs of subsystems are not 
shown.). States that belong to two of these sets are at most 
bipartite entangled, and states that belong to all three sets are 
separable. Any convex sum of bipartite entangled states (de- 
picted by T2) is considered at most bipartite entangled, even 
though the state might not be separable with respect to any 
bipartition. Similarly, any convex sum of tripartite entangled 
states (depicted by rs) is considered at most tripartite entan- 
gled. Only states that can not be obtained as a convex sum of 
at most tripartite entangled states are four-partite entangled. 

fc-partite entanglement in mixed states. We first start out 
describing the underlying idea in rather general terms, 
followed by a specific realization tfiat satisfies all of the 
desired properties. What we aim at is a set of functions 

n/2 

rfc,„(g) = /(£<)-^af'")^/„-(^,) (2) 

i=l j 

defined in terms of functions / and fij where the index 
j labels all inequivalent bipartitions [18] of the n-partite 
system in an i-partite and an (n — i)-partite component, 
referred to as i-bipartitions in the following. Further- 
more, the functions / and fij, and the scalar weight fac- 
tors flj-*''"'' need to satisfy the following conditions 

I Tk^nig) is convex, i.e. J^iPi'^iod > '^(Z^jKft) 
II fig) > and /„(g) > 0, V^, j, g>0. 

III /(?/;) = /ij(V') if IV-") is bi-separable with respect to 
the j-th i-bipartition. 

IV af'"^ > 0, yi,k,n. 

Convexity of tj, „ allows us to restrict the following dis- 
cussion to pure states: if Tk^n is non-positive for all pure 
states with less than fc-partite entanglement, condition I 
entails that a positive value of Tk^n identifies (at least) 
fc-partite entanglement in mixed states. What remains 
to be done is to tailor the prefactors a^'"^ in order for 



Tk^n to have the desired properties for pure states. Since 
fiji"^) coincides with /(vE*) if l^*) is separable with re- 
spect to the j-th i-bipartition, it is sufficient to charac- 
terize the separability properties of pure fc-partite entan- 

(k n) 

gled n-partite states, and choose the weights a] ' such 
that I]ij|f,s '"''/»j (*fe'n) > /(*fc'n) for any fc'-partite 
entangled state with k' < k, where the sum runs over 
all bipartitions with respect to which l^'fc'n) is separable. 

Due to the positivity of fij and a^'''"' this directly im- 
plies that Tfc „ is non-positive for all states which are not 
at least fc-partite entangled 

As indicated in Fig. 2, pure states with only a small en- 
tangled component are biseparable with respect to many 
bipartitions, so that many components fij{'9) coincide 
with /(^) and the weights can be chosen comparatively 

(k 71) 

small. Choosing the weight factors ' increasing with 
fc will thus allow us to arrive at the desired hierarchies. 

At the bottom of the hierarchies lies r2,„ which has to 
be non-positive for all completely separable states I'ifi^n)- 
Since I'i^i.n) is separable with respect to any biparti- 
tion, we have /(^'i,™) = fij{^i,n) Vi,j due to (II). Any 
choice satisfying il^'"'' — 1 will therefore result in 
T2,n{'^i,n) = for any completely separable state |^i,n). 
In order to proceed we need to tailor the a^'"^ such 
that T3^„ is non-positive for all pure states that contain 
less than tripartite entanglement. As depicted in Fig. 2 
with the exemplary case of n = 5, any pure state |5'2,n) 
(for n 7^ 4) is separable with respect to at least m 2- 
bipartitions, where m is the largest integer < n/2 [19]. 
Accordingly, a^'"' — Si2/'m is a valid choice for T3,„. 

Typically there is not a unique choice for the weights 
^(fe:")^ and the resulting freedom can be used to opti- 
mize the functions ti^^ ^ 

for specific quantum states. As 
a rough rule of thumb we found that for states with 
highly mixed reduced density matrices choices with large 
weight factors a,-'^'"-' for i ~ n/2 and small or vanish- 
ing ones for i <^ n/2 yield strong criteria. For example 
for odd n > 5, a^'"'' = Sis/m and ap'"'' = ^ii/m are 
both valid choices to define T3_„, but we found the for- 
mer to result in a stronger criterion. Similarly, for n > 8 
^(3,n) _ ^^^/(jji(j(i _ l)/2) typically leads to an even 
stronger criterion. Since picking a good choice for the 
weight factors helps to identify good criteria, we refrain 
from providing a systematic description for the construc- 
tion of the ai''^^, but rather depict choices for n < 7 that 
we found to yield good results in Table I. The functions 
'Tk,n with these specific coefficients then define a full hi- 
erarchy of necessary separability criteria for any system 
size n. 

As is the case for any attempt to detect entanglement 
beyond 2 x 3-dimensional systems [6], a tool can either 
identify entanglement or it can identify separability, but 
there is none that can assert with certainty whether a 
state is entangled or separable. Also here a non-positive 
value of Tfc^n does not necessarily imply that the con- 
sidered state was not fc-partite entangled, but it could 
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K = 2 


[0,0,0,0,0,1/462] 


[0,0,0,0,1/462] 


[0,0,0,0, 1/126] 


[0,0,0,1/126] 


[0,0,0,1/35] 


[0,0,1/35] 


[0,0,1/10] 




K = 3 


[0,0,0,0,0,1/10] 


[0,0,0,0,1/10] 


[0,0,0,1/10,0] 


[0,0,0,1/6] 


[0,0,0,1/3] 


[0,0,1/3] 


[0,1/3,0] 




K = 4 


[0,0,0,0,0,1/3] 


[0,0,0,0,1/3] 


[0,0,0,1/3,0] 


[0,0,1/3,0] 


[0,0,1/2,1/3] 


[0,0,1/2] 


[0,1/3,1] 




K = 5 


[0,0,0,1/3,0,1/3] 


[0,0,0,1/2,1/3] 


[0,0,0,1/2,0] 


[0,0,1/3,1/2] 


[0,0,1/2,1] 


[0,0,1] 


[0,1,1] 




K = 6 


[0,0,0,1/3,1/2,1/3] 


[0,0,0,1/2,1/2] 


[0,0,0,1/2,1] 


[0,0,1/3,1] 


[0,0,1,1] 


[0,1,1] 


[1,1,1] 




K = 7 


[0,0,0,1/3,1/2,1] 


[0,0,0,1/2,1] 


[0,0,0,1,1] 


[0,0,1,1] 


[0,1,1,1] 


[1,1,1] 






K = S 


[0,0,0,1/3,1,1] 


[0,0,0,1,1] 


[0,0,1,1,1] 


[0,1,1,1] 


[1,1,1,1] 




[1,1] 


K = 5 


K = 9 


[0,0,0,1,1,1] 


[0,0,1,1,1] 


[0,1,1,1,1] 


[1,1,1,1] 




[1,1] 


[0,1] 


K = 1 


K = 10 


[0,0,1,1,1,1] 


[0,1,1,1,1] 


[1,1,1,1,1] 




[1] 


[0,1] 


[0,1/2] 


K = 3 


K = 11 


[0,1,1,1,1,1] 


[1,1,1,1,1] 




[1] 


[1/3] 


[0,1/3] 


[0,1/10] 


K = 2 


K = 12 


[1,1,1,1,1,1] 






N = 2 


iV = 3 


iV = 4 


N = 5 





TABLE I. Specific choices for the weight factors a\ that define valid criteria tkn to detect /f-body entanglement in an 
A'"-body system. Each vector in the table contains the elements [ffl^^', ■ • • , a^^^'] (for even A'") and [a!^^^\y2^ • ■ • , ii^^''] (for 
odd A''). The upper left half ranges from A = 12 to A^ = 6, the lower right from A = 2 to A = 5. 




FIG. 2. Schematic representation of two pure five-partite 
states: a bipartite entangled (top) and a tripartite entangled 
(bottom). The bipartite entangled state is separable with 
respect to the 2-bipartitions S1S2 — 82.8485 and 818285 — 
<S3<S4. Indeed, any pure bipartite entangled five-partite state 
is biseparable with respect to at least two 2-bipartitions. Pure 
tripartite entangled five-partite states on the other hand can 
be separable with respect to one 2-bipartition only. 



also be due to the fact that rfc^„ is not strong enough to 
identify the targeted entanglement in the specific state. 
If the latter is the case, one can improve Tk,n provided 
there are additional properties that can be exploited. 
Here we would like to demonstrate this with the exam- 
ple of VF-states [20], i.e. states with a single excitation 
\W) = X^i'^^d*), where \i) is a short hand notation for 
the state with the i-th subsystem in its excited state and 
all other subsystems in their ground state. These states 
attract particular attention since they occur naturally 
in excitation transport processes [16, 21], and they also 
permitted the observation of genuine multi-partite entan- 
glement of an eight ion string [22]. 

If such a ly-state is bi-separable with respect to an 
i-bipartition, so that \W) = |0i) ig) |02), then one of 
the components \4>i/2) needs to be completely separa- 
ble, because otherwise there would be a finite ampli- 
tude for two excitations [23]. Consequently, bisepara- 
bility with respect to an i-bipartition {i < n/2) implies 
biseparability with respect to at least i one-bipartitions, 
i(i — l)/2 two-bipartitions, and similarly for larger bi- 
partitions. These additional separability properties per- 



mits to identify significantly lower values for the weights 

(k n) 

al ' than those for general states as given in Table I. In 
contrast to above, where we found strong criteria based 
on bipartions of ~ n/2 subsystems, in the case of In- 
states it is rather advantageous to focus on 1-bipartitions: 
the weights a-'^'"-' — Sii/{n — {k — 1)) for A: ^ 2 and 

(k n) 

a-i ' — dii/n for k — 2 provide a strong hierarchy 
for VF-states. In a similar fashion, the present criteria 
can also be adjusted for different classes of states, such 
as more general Dicke states [24] or potentially states 
with permutation symmetries [25] . 

So far we have discussed the hierarchies in a rather ab- 
stract setting, assuming the existence of functions that 
satisfy the above list of properties I to IV. Let us be- 
come more specific now and present a possible choice of 
such functions. It is based on the fact that a twofold 
tensor-product ® [^f) of a state with itself features 
very specific invariance properties if l^*) is not genuinely 
n-partite entangled [26]: |^') is biseparable with re- 
spect to a biseparation that divides the system in the 
components A and B if and only if the two-fold state 
l^*) (g) l^f) is invariant under the permutation that per- 
mutes the two j4-components (or, analogously, the B- 
components). Taking / to be a function g of l^*) (E) \'^) 
and fij = g{Ilij\'^) (g) l^*)) where Hij is the permutation 
that permutes the A-components associated with the j- 
th i-bipartition makes sure that condition III is satisfied. 
Condition I, i.e. convexity of r^^n, is in general difficult 
to achieve, but 

/ = VWln^^M^ and 
= ^J{^s\Tiij0(SgliU^s) (3) 

with the global permutation 11 and a product vector 
\^s}: are convex resp. concave [8, 9]. 

For pure states / coincides with ^ {^s\g'E q\^s) 
d'l') (g) 1^) is invariant under 11), so that condition III is 
satisfied, and the present specific choices for / and fij are 

indeed non-negative. As long as a^'^'"'' are non-negative 
as they should be according to condition IV, Eqs. (2) and 
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(3) with the weight factors al such as those given in 
Table I provide a vahd reahzation of a hierarchy following 
conditions I through IV. 

We have tested the performance of these hierarchies for 
different, exemplary cases, comparing it with previously 
known criteria [4, 27] for 4- and 6-partite spin-squeezed 
states and 4-partite W-states. This comparison is shown 
in the section A of the appendix. This test demonstrates 
how the hierarchies Tk^n and t™„ often outperform prior 
techniques, especially in presence of strong mixing. This 
is remarkable, in particular, since existing criteria have 
been specifically tailored to address entanglement prop- 
erties of a given class of states, whereas we just varied the 
coefficients a|'^'"'^ retaining the same analytic form. By 
a systematic application of the hierarchy Tk^n to a given 
system of interest it is possible to gain a deep insight in 
its entanglement structure, as noticeable for the exem- 
plary case of spin-squeezed states in the appendix, where 
a rich underlying multi-partite entanglement landscape 
is uncovered in a broad interval of the spin-squeezing 
parameter where previously only bipartite entanglement 
was detected. 

As argued above, the possibility to detect /c-partite 
entanglement for varying k also provides a refined in- 
sight in entanglement dynamics. This is substantiated in 
Sec B of the appendix with the investigation of a fully 
connected 12-partite graph state undergoing a dephasing 
evolution, where it is emphasized how the different types 
of fc-partite entanglement decay on different time scales. 

The possibility to explore the dynamics of arbitrary fc- 
partite entanglement in turn enables to uncover relevant 
physical features of a given system, as we exemplify with 
the verification of three-body interactions [28] in section 
C of the appendix. 

These examples underline the usefulness of the present 
hierarchies; the specific framework of permutation oper- 
ators that we have used for the explicit construction of 
the present hierarchy is by no means the sole way to ar- 
rive at such a hierarchy. Many other typically employed 
tools, such as entanglement witnesses [29] or positive 
maps [30] bear potential for a systematic construction. 
Also, whereas we have focussed here on the classification 
of fc-partite entanglement, a classification in more refined 
classes according to LOCC-inequivalence [20, 31] seems 
feasible. 

Inspiring discussions with Lukasz Rudnicki, Otfried 
Giihne and generous financial support by the European 
Research Council are gratefully acknowledged. 



I. APPENDIX 

This appendix shows various applications of the hierar- 
chies Tk^n in order to provide a quantitative estimation of 
their performance and to exemplify their potential uses. 
In the first section we compare the strength of the hi- 
erarchies with previously introduced separability criteria 
[4, 27] for W-states in I A I and spin-squeezed states in 



I A 2. Subsequently, in section IB, we utilise one hier- 
archy to investigate the decay of multi-partite entangle- 
ment for the case of a genuine 12-partite entangled state 
undergoing a dephasing dynamics. Furthermore in sec- 
tion I C we demonstrate how detecting fc-partite entan- 
glement in an n-partite system with k ^ n provides a tool 
to distinguish between three and two body interactions. 

A. Comparison with previous criteria 

1. W-states 

A criterion to detect fc-partite entanglement in 4- 
partite M^-states, i.e. states with a single excitation, has 
been described in [27] based on expectation values of a 
given density matrix g with respect to the four states 

\Wi) = ^(11000) +e'''^i 10100) +e*'^^ 10010) -he*"^^ 10001)) , 

\W2) = ;^(|1000) -e^'^HOlOO) -e*"^^ 10010) -he*"^^ 10001)) , 
\ (4) 
|W^3) = -(11000) +e''^i 10100) -e*"^^ 10010) -e*-^^ 10001)) , 

\Wi) = ^(11000) -e* 10100) +e*'^^ 10010) -e*-^^ 10001)) . 

As shown in [27] , the variance based quantity 

4 

l^{g)^l-Y,{WMW.)\ (5) 

i=l 

permits to identify bi-, tri- and four-partite entangle- 
ment: if A(gi) < Afe for some choice of the phases (^i, 
(/>2 and 03, with A2=3/4, A3=l/2 and A4=5/12, then g 
is fc-partite entangled. 

In order to compare the entanglement criteria we 
generated ensembles of random W-states g{A) = 
AA^ /Tt:{AA'^) in terms of random matrices A whose ele- 
ments are independently and normally distributed. This 
ensures a distribution according to the Hilbert-Schmidt 
measure [32] , and the distribution of the entropies of the 
states p{A) is determined by the ratio of the mean p, and 
variance of the Gaussian distribution from which the 
matrix elements are obtained. With states generated in 
this fashion we will compare in the following 

dk{g) = min A(gi) - Afc, (6) 

i.e. the criterion from [27] with an optimal choice 
of phases with the hierarchy t^„, where we also 
optimise over the separable states \^s) introduced in 
Eq.(efeq:tau). 

The choices /i = and cr^ = 1 lead to an ensemble of 
intermediate mixing with an entropy distribution shown 
in blue in the inset of Fig. 3. All states are detected 
as bipartite entangled according to both ^2 and r^4, re- 
gardless of their mixedness. For the case of tripartite 
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entanglement, however, there is a rather striking differ- 
ence in the performance of the two criteria as it can be 
seen in Fig. 3, where is plotted against t^^ for this 
ensemble. Whereas detects only 5% of the ensemble 
members as tripartite entangled, t^^^ manages to detect 
65%, and there is no state detected by ds that is not 
detected by t^^. 

In the case of genuine 4-partitc entanglement it had 
already been observed that t^4, i.e. the genuine mul- 
tipartite entanglement criteria presented in [33], is not 
particularly strong for W^-states, and, indeed, there are 
various cases where d^ performs better than t^^: 0.05% 
states yield a negative value of c?4 whereas only 0.003% 
lead to positive t^^^. This can however be accounted for 
by a modification of (Eq.(III) in [33]) that is better 
suited for ly-states. With this modification the topmost 
members of both hierarchies detect 0.05% of all states 
and agree in more than 90% of the detected cases. 

Inspecting the entropies of the states that are detected 
by r;^4 but not by ^3, one realizes that performs bet- 
ter in particular for rather highly mixed states. We there- 
fore considered two additional ensembles: with jl — 3 (20) 
and = 2 (8) what results in rather high (low) entropies 
as depicted in red (yellow) [34] in the inset in Fig. 3. For 
close to pure states (yellow distribution) both criteria 
perform equally well, but in the case of highly mixed 
states (red distribution), the situation changes drasti- 
cally: (^3 does not detect any state, whereas r™4 still 
manages to identify 5% and thus proves to be signifi- 
cantly stronger for highly mixed states. 

To better explore the different performances of the two 
criteria depending on the mixing of the quantum states 
we consider now a Werner-type state 



Qq = q\Wi){Wi\ + {l-q) 



2" ' 



(7) 




FIG. 3. Comparison between the hierarchy and bound 
Eq.(5) for the case fe=3. In the inset the distribution of en- 
tropies for the three different batches of data are depicted. 
The main picture shows obtained values for the two criteria 
with states whose entropy follows the blue (middle) distri- 
bution in the inset, each point representing a single random 
mixed state. All states under the 1/ = axis violate bound 
Eq.(5) and are therefore detected by d^. Tripartite entan- 
glement is detected by the presented hierarchy when t^4 is 
positive. No state with d3<0 leads to a non-positive value of 
T^4. The majority (~60%) of random states, especially for 
higher mixedness, yield a positive r^^i, but a variance A(p) 
not small enough to detect tripartite entanglement according 
to Eq.(5). The value for rf^^ has been normalised to unity for 
a perfect W state. 



following inequality 



Var(J,) > -F, 



fc/2 



U/2 



(8) 



in order to systematically explore mixedness via a single 
parameter q. With both inequality Eq.(5) and the hier- 
archy r™„, 4-partite entanglement is detected for q>2/3, 
and bipartite entanglement for all values of q>Q. How- 
ever, the threshold value for tripartite entanglement de- 
tection is found to be g « 0.58 by violation of Eq.(5), and 
(7=0.3 through the hierarchy. That is, r^4 detects nearly where 
twice the parameter range to be tripartite entangled than 



where A; is a divisor of n, and Jz are collective spin 
operators and is a set of numerical bounds that can be 
found in [4]. Violation of Eq.(8) constitutes a sufficient 
criterion for k + 1-partite entanglement, and we compare 
its performance with the hierarchy of criteria Tfc_„ given 

by the choice of coefficients a^'^'"'' listed in Table I. 



Spin squeezed states are of the form 



Ix) 



$0) 



V2 



(|o.) + IU) 



(9) 



(10) 



2. Spin-squeezed states 

Also for spin-squeezed states [35] a criterion for the 
identification of fc-partite entanglement has been de- 
scribed before. Spin-squeezed states are entangled, and 
the strength of spin-squeezing is related the degree of en- 
tanglement [4] . More specifically, the collective spin of a 
composite system comprised of n spin-1/2 particles which 
is fc-partite entangled (fc=l means separable) satisfies the 



is the eigenstate of Jx to the eigenvalue n/2. Since we 
are interested in mixed states, we consider the admixture 
of white noise, i.e. Werner-type states: 



Qx,r, = V\x){x\ + (1 -^7) 



1. 



(11) 



with the n-partite identity 1„. For the present com- 
parison we consider the cases n — 4(6) and gener- 
ated 4000 (2500) random states each with x and 77 
distributed uniformly in the intervals 0<x<0.5 and 
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FIG. 4. Comparison of the spin-squeezing inequality Eq.(8) and the hierarchy r^^^ for the case of n=4 (left) and n=6 (right) 
spin 1/2 particles. Each point corresponds to a single random state. The solid curves represent the bound Eq.(8) for k=l,2 
and 4 (n=4) and 1,2,3 and 6 {n—6) respectively from top to bottom. Lying below a certain k curve gives a sufficient condition 
for k + 1-partite entanglement. Colour-code: states which lead to a rfc_„ > are brown (fc=6), cyan (fc=5), magenta (fc=4), 
green (fc=3), red (fc=2); blue points correspond to states for which no entanglement has been detected. The curve for fe=3 in 
n=6 cannot be rigorously generated for all values of (Jx)' the solid line depicts where Eq.(8) is saturated, and the dashed line 
its analytical continuation which overestimates the bound. 



0.8<ry<l. To achieve optimal performance of the spin 
squeezing inequality Eq. (8), expectation values of jTz 
are not taken vifith respect to Px,»7' with respect to 
exp^—iiyj^x) Qx,ri exp(iz^j7a;), and an optimisation over the 
angle v is performed. 

Fig. 4 depicts the minimal variance Y&t{J^z) as function 
of J^z for various values of k with n — 4, and n = 6. The 
black lines delimit the areas in which Eq. (8) is violated, 
and A: -|- 1-partite entanglement is identified; in some cases 
different degrees are not resolved due to the limitation 
that fc be a divisor of n. Each of the colored points 
corresponds to one random state, and the color coding 
refers to the degree of entanglement that is detected by 
the positivity of Tk^n (details in the caption of Fig. 4). 

For n = A, the hierarchy rfc,„ is strictly better than 
Eq. (8): not only is it able to differentiate between 
tri- and four-partite entanglement, but within the set of 
states detected as bi(tri)partite entangled by the spin- 
squeezing inequality it shows that more than half of 
them (^ 60% (52%)) are actually /c-partite entangled 
with k > 2(3). Furthermore, there are no states for which 
Tfe^4 detects a smaller number of entangled components 
than as estimated via the spin-squeezing inequality. 

Since for n = 6 there are a few states where the vio- 
lation of Eq. (8) leads to an estimated /c-partite entan- 
glement larger than detected by the present hierarchy, 
we shall analyze this case in more detail. Out of all the 
states detected as bipartite entangled by T2_6j and de- 
picted in red in Fig. 4, 40% are non-squeezed (they lie 
above the k=\ curve) and, therefore, can not be detected 
as entangled by Eq. (8) as a matter of principle. 53% of 
the states depicted in red are detected as bipartite en- 
tangled also via their spin squeezing properties, but for 
7% Eq. (8) detects a larger degree of entanglement than 
Tk^n as depicted by the red dots below the fc=2,3 curves. 
Out of the states detected as tripartite entangled by the 



present hierarchy (depicted in green in Fig. 4), 66% lie 
between the k=l and the fc=2 curve and are detected 
only as bipartite entangled by Eq. (8); ^ 18%, however, 
are detected as four-partite entangled, since they lie be- 
low the fc=3 curve. 

About 52% (49% considering only squeezed states) of 
the entire sample is detected as 4-partite, 5-partite or 
6-partite entangled by Tk^n (depicted by magenta, cyan 
and brown respectively in Fig. 4), whereas Eq. (8) does 
not permit to distinguish between 4-partite entanglement 
and fc-partite entanglement with fc > 4. In addition, a 
substantial portion of those states is detected by Eq. (8) 
as bipartite entangled only, or is not detected as entan- 
gled at all. 

After all, Eq. (8) and the present hierarchies result in 
inequivalent assessments, and the result of a compari- 
son likely depends on the sample of states chosen. Here 
we have been choosing the parameters such that there is 
a substantial portion of states with significant squeezing, 
and found that Eq. (8) identified a larger degree of entan- 
glement in ^ 7% of the squeezed states, whereas among 
^ 47% of those states r/j^g detects a larger number k 
of entangled components. That is, despite the fact that 
Eq. (8) has been designed specifically for spin-squeezed 
states, and un-squeezed states have not been taken into 
account in the comparison, our present hierarchies out- 
perform spin-squeezing inequalities significantly more of- 
ten than vice versa. 



B. Decoherence of open quantum systems 

So far in this appendix the focus has been on compar- 
ing the detection strength of hierarchy Tk^n with previ- 
ously known criteria for specific classes of states. Now we 
aim to show the potential of the hierarchies to identify 
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physical properties of multi-partite entanglement. This 
is exemplified by monitoring the qualitative changes of 
entanglement induced by decoherence in a multi-partite 
system undergoing a dephasing evolution. 

To this end we choose a system of n = 12 two-level 
systems and consider the impact of dephasing on a fully- 
connected graph state [36] , where each subsystem is sub- 
ject to a dephasing channel with Kraus operators Ei = 
V(l+exp(-7t))/2 1 and E2 = ^(1 - exp(-7t))/2 a„ 
with 7 the single particle decoherence rate. With these 
single-particle channels, the multi-partite dynamics reads 



g{t)= J2 



E, 



(12) 



TKN 




Again an optimisation over the product vectors |<1>5) in- 
troduced in Eq.(3) of the manuscript is performed in or- 
der to detect entanglement as reliably as possible with the 
hierarchy Tk^n and the coefficients a^'^'"'' from Table I. In 
the present case we verified for n < 8 that 



(|0)-|-i|l))^"®(|0)-i|l))^"/2" 



(13) 



maximizes Tfc^„ for all times independently of k and n, 
and consequently used this choice of \^s) also for n > 8. 
The decay properties depicted in Fig. 5 are thus known 
as analytic functions of t. 

Fig. 5 shows the decay of rfe^„ for n = 12 and k ranging 
from 2 to 12, where the substantially different behavior of 
fc-partite entanglement for different k is noticeable. The 
instant tg at which Tfc^„ becomes negative, so that no k- 
partite entanglement is identified anymore, is depicted 
in the inset as a function of k. As can be seen in this 
specific example 12-partite entanglement can be observed 
until ~ 0.12 only, but the life-time of 10-partite entan- 
glement is already more than 2.5 times longer. Further- 
more, the lifetime of bipartite entanglement (7^ ~ 2.61) 
exceeds that of genuine 12-partite entanglement by a fac- 
tor of ~ 22, i.e. by more than an order of magnitude. 
The very different dynamical behavior that the various 
fc-partite entanglement feature shows that few-partite en- 
tanglement can behave very differently from multi-partite 
entanglement. 



FIG. 5. Entanglement decay for an n = 12-partite system 
initially prepared in a fully-connected graph state and evolv- 
ing according to a dephasing channel for k ranging from 2 to 
12 depicted with increasingly long dashing. The inset shows 
the instants ts from which no fc-partite entanglement is de- 
tected anymore. 
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FIG. 6. Comparison between the exponential growth of 
'''3,o{\tp{'t)}) in the case of a two-body (red) or three-body 
(blue) nearest neighbour spin-spin interactions as given in 
eq.(14). Solid lines correspond to open boundary conditions, 
i.e. a chain; dashed lines to periodic boundary conditions, i.e. 
a ring. For short times the growth is monomial in time cx 
and the exponent is clearly shown to depend only on the type 
of the interaction, with all two-body interaction converging 
to Q=2 whereas a=l for three-body interactions, such that a 
allows to discriminate between them. 



C. Identification of three-body interactions 



Finally, as an example of how detecting fc-partitc en- 
tanglement can provide insight on the dynamical features 
of a given system, let us demonstrate how the present 
hierarchies can be used to verify the existence of three- 
body interactions whose engineering is currently actively 
debated [28, 37]. 

Since typically more pairwise interactions than inter- 
actions among triples are required to create a /c-partite 
entangled state, the onset of entanglement growth start- 
ing from a separable state can be expected to be a dis- 
criminator between these two types of interaction. For 



sufficiently short times entanglement grows in a mono- 
mial fashion Tk.n{\'i'{t))) (x t" and, for a given system 
size n and a suitably chosen fc, the exponent a is a signa- 
ture that discriminates two and three-body interactions, 
independently of their strength or specific form. This is 
shown here with the exemplary case of tripartite entan- 
glement in a five-body system (fc=3, n=5). 

We consider the dynamics generated by Hamiltonians 
of the following form: 
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(14) 

k 

where Ui, aj, cr£ are taken randomly from the three Pauh 
matrices and Afc are random couphng coefficients. A 
specific reahzation of these Hamiltonians will then gen- 
erate the time evolution of an initially separable state 
|V(0))- 111111). 

Fig. 6 shows how indeed the growth of T3.5 is monomial 
in time, as the time derivative of Logio(T3,5) approaches 
a constant value for t — )■ 0. The curves corresponding to 
two-body interaction Hamiltonians (i?26) are depicted in 
red, those originating from three-body interactions (H^h) 



in blue: it is neatly shown how the exponents a converge 
to a value which depends exclusively on the nature of the 
interaction, independently of its details. For all two-body 
Hamiltonians a converges to the value a— 2, whereas for 
all the three-body interactions a converges to 1. Only 
for larger times signatures of the specific realisation be- 
come apparent, as deviations from the asymptotic val- 
ues. This approach therefore provides a reliable method 
to witness whether a given interaction is of two or three- 
body nature. What is shown in Fig. 6 is in sharp con- 
trast to bipartite entanglement, which grows linearly in 
time both for two-body and three-body interactions and 
is thus of no help. This emphasises how different tasks 
may require probing fc-partite entanglement for 2<k<n, 
and that the access to various fc-partite entanglement al- 
lows us to grasp underlying dynamical features. 
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